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Abstract 

The instantonic approach for a </> 4 model potential is reexamined in the asymptotic limit. The 
path integral of the system is derived in semiclassical approximation expanding the action around 
the classical background. It is shown that the singularity in the path integral, arising from the 
zero mode in the quantum fluctuation spectrum, can be tackled only assuming a finite (although 
large) system size. On the other hand the standard instantonic method assumes the (anti)kink 
as classical background. But the (anti)kink is the solution of the Euler-Lagrange equation for 
the infinite size system. This formal contradiction can be consistently overcome by the finite size 
instantonic theory based on the Jacobi elliptic functions formalism. In terms of the latter I derive 
in detail the classical background which solves the finite size Euler-Lagrange equation and obtain 
the general path integral in finite size. Both problem and solution offer an instructive example of 
fruitful interaction between mathematics and physics. 
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I. Introduction 



The tunneling through a potential barrier plays a central role in many areas of the 
physical sciences. It is well known that quantum tunneling is an intrinsically nonlinear 
phenomenon which cannot be attacked by standard perturbative methods. As an example, 
take a symmetric double well potential (DWP) often encountered in the form of a two 

nnn 

level system in condensed matter physics [1|, I2J, |3|]: the stable minima are asymptotically 
connected by the classical (anti)kink solution, a charge conserving domain wall, whose energy 
is inversely proportional to the quartic force constant (5 > 0) of the 4 potential at a given 
vibrational frequency inside the well. Then, even small quartic nonlinearities induce large 
classical energies which, in turn, enhance the effects of the quantum fluctuations in the 
overall probability amplitude for a particle to move from one minimum to the next. This 
explains why standard perturbative techniques fail to describe tunneling processes while 



semiclassical methods are known to provide an adequate conceptual framework 

In general, 4 models also host potentials having a negative sign of the quadratic term 
(5 < 0) [6]. In this case the mid-point of the potential valley would be a classically and 
locally stable configuration, a metastable state from which the particle would tend to escape 
via quantum tunneling. In such state, known in quantum field theory as "false vacuum" [7|, 
the system has a finite lifetime due to the presence of a negative eigenvalue in the quantum 
fluctuation spectrum which governs the decay into the "true vacuum". The classical solution 
of the Euler-Lagrange equation for metastable potentials is a time-reversal invariant bounce 



is 



carrying zero charge whose energy is inversely proportional to the anharmonic parameter 
5. Then, also analysis of metastable states and decay rates in tunneling problems are 
unsuitable for perturbative methods. 

Semiclassical methods usually extend the WKB approximation lOj and permit to evalu- 
ate the spectrum of quadratic quantum fluctuations around the classical background. While 
such methods had been first applied to infinite size systems, extensions to finite size sys- 
tems have been later on developed H[ to obtain correlation functions and spectra of (/> 4 
and other model potentials [ijj]. Classical mechanical systems perturbed by 

spatio-temporal noise are known to be sensitive to finite size effects and, in particular, there 
is growing interest in metallic nanowires and micromagnets with magnetization reversal, 
in view of their technological potential. Intense research is in progress also on biological 
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systems and specifically DNA to establish the finite size effects on the physical properties 

m HQ- 

While a detailed study of decay rates in finite size metastable potentials has been pre- 
sented in Ref. jl9|. this paper reconsiders the infinite size instantonic theory for a bistable 
4> A potential: in particular, I will point out the feature which justifies a consistent finite size 
analysis entirely based on the Jacobi elliptic functions formalism. The latter is not exactly 
well known among the physicists. A stimulating introduction to the elliptic functions has 



been published by Good [20(] while a systematic treatment can be found in the Refs. |2l|, |22|. 

Thus, I feel, this is an interesting situation in which physics and mathematics mutually 
benefit by each other: the well known instanton problem offers the ground for an application 
of the Jacobian functions while the latter provides the fundamental tool to extend the 
instanton theory to cases of practical interest. 

The key point of the instanton problem can be synthesized as follows: due to the time 
translational invariance a Goldstone mode appears in the system. As shown in detail below, 
such mode is proportional to the (anti)kink velocity and represents the ground state for the 
quantum fluctuation spectrum. Accordingly, the associated zero mode eigenvalue Eq breaks 
the Gaussian approximation and makes the Euclidean path integral divergent. The trouble 
is overcome by regularizing the fluctuations determinant and evaluating the contribution of 
the zero mode separately. In doing that however, the size L of the system (the period of 

— 1/2 

the instanton along the time axis) is assumed large but finite so that e oc L and the 
path integral is finite. Then, to be consistent, also the classical equation of motion should 
be solved assuming a finite L and the (anti)kink, peculiar of the infinite size theory, should 
be recovered as a limit case of the more general, finite size (anti) instanton solution. This 
problem is discussed in detail in Section II. In Section III, I revisit the standard Euclidean 



path integral 23|, |24( implemented by the functional determinants method which permits to 



derive, in the asymptotic limit, the contribution stemming from the quantum fluctuations 
around the (anti)kink 25, [26]. A consistent extension of the functional determinants theory 
to the finite size system allows one to obtain the path integral for the finite size instanton. 
The general path integral formula are given in Section IV while Section V contains some 
final remarks. 
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II. Classical Equation of Motion 



Take a DWP described by: 



V(x) = —x 2 + -x 4 + -a 4 (1) 

with minima located at x = ±a (V(x = ±a) = 0) and positive 5 in units eVA~ 4 . Note 
that the relation a = \jMu 2 J8 holds thus posing a constraint on the potential parameters. 
Eq. (pQ) is equivalent to: 



= ^-(^-a) (ar + o) . (2) 

At a fixed a, a larger oscillator frequency uj inside the well leads to an increased non- 
linearity. Obviously there cannot be classical motion through the barrier of a DWP and, 
in the real time formalism, the very solutions are given by the vacua: x c i = ±a. However, 
applying a Wick rotation which maps the real time t onto the imaginary time r, the DWP is 
turned upside down and the classical equation of motion admits a non trivial solution j^ . 
The latter is the classical path x c i(t) which minimizes the action in the imaginary time path 
integral formalism. Let's find it solving the Euler-Lagrange equation which reads: 



Mx d (r) = V\x d {r)) (3) 
where V means derivative with respect to x c \. Integrating Eq. Q, one gets 



^x 2 d {r)-V{x cl {r))= E (4) 

with integration constant E. Physically, we are looking for a path covering the distance 
between —a and +a in the period L. This means that the solution has to fulfill the boundary 
conditions: 



x d (r = ±L/2) = ±a. (5) 
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Moreover, from Eq. (j4j) it is easily seen that the relation 



x d (r = ±L/2) = ^2E/M (6) 

holds on the r— range boundaries. Hence, only (anti)instantons associated to values 
E > can describe the physics of the DWP. This marks an essential difference with respect 
to the bounce solution of the Euler-Lagrange equation for metastable potentials which, in 
fact, exists only for E < Of course there is a deep physical reason behind such 

difference in sign: the instantons fulfills antiperiodic boundary conditions and carries finite 
topological charge whereas the bounce is a time reversal invariant object fulfilling periodic 
boundary conditions. 

Integrating Eq. one gets: 



[M r a(T) dx 

T-T = ± 




2 J Xcl(T0) ^E + V{x) 

(7) 

To represents the time at which the middle point of the valley in the reversed potential is 
crossed. As the system is translationally invariant, tq can be placed everywhere along the 
imaginary time axis, r G [— L/2, L/2]. This arbitrariness ultimately represents the physical 
origin of the zero mode in the stability equation which governs the quantum fluctuations 
spectrum (see below). Eq. ([7j) admits a solution which interpolates between the vacua 
providing evidence for the existence of quantum tunneling. Such solution is in fact a family 
of paths whose shape depends on the associated energy E. 

To be specific, it is clear from Eq. (J7J) that the zero energy configuration (E = 0) is 
consistent with a physical picture in which the particle starts at x(t = — oo) = —a with 
zero velocity, falls in the valley and reaches the top of the adjacent hill at x(t = +oo) = +a 
with zero velocity. Integration of Eq. (|7j) with E = leads to the classical path: 



x c i{ T ) = ±atanh[o;(r — Tq)/v2]. 



(8) 



This is the well known (anti)kink solution which exists if and only if the potential minima 



are asymptotically connected. In the spirit of the Matsubara formalism 



24 



281 ]. the period 



L can be mapped onto the temperature axis as L = H/KbT. Hence, for practical purposes, 
Eq. (jSJ) holds in the low T / large L range whose upper / lower bound is set by r* such that 
\ x ci{ T *)\/ a — 1- As long as the latter condition is fulfilled the (anti)kink solution may be 
assumed as a reliable approximation for the finite size (anti)instanton. This terminology, 
born in quantum field theory {29! . I30] ]. stems from the fact that the tunneling transition 
happens in a short time, almost instantaneously. From Eq. ([8]) one also notes that potentials 
with larger 10 sustain classical kinks in a broader temperature range. 
Take now E > and define: 



/ \ 1 X cl{r) 

Xd[T) = T~ 

K =2 + «' < 9 » 



Then Eq. ([7]) transforms into: 



1 

T - T = ±- 



dx 



(10) 



'Xci (to) Vx* ~X 2 + «/2 
which requires some work to be solved in terms of Jacobian elliptic functions. The integral 
in the r.h.s. of Eq. (fIU|) yields: 



dx 



VX 4 ~X 2 + k/2 



/ 4 2 
6= 2 arctan ( \ —Xd 

2 V V2^ 



-F(9,s) 



(11) 



F(8, s) is the elliptic integral of the first kind with amplitude 9 and modulus s: 



m*)= [ ~r ==, (12) 

whose quarter-period is given by the first complete elliptic integral K(s) = F(jr/2, s) 21]. 
In terms of F(6, s), one defines the Jacobi elliptic functions sn-amplitude, cn-amplitude, dn- 
amplitude : 
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sn[F(6,s)] = sm9 
cn[F(6,s)] = cosd 



dn[F(9,s)] = Vl-s 2 sin 2 e. (13) 
From Eqs. ffTUl) . (ITTj) one gets: 



±\I-2u(t-t )= F(e,s) (14) 



and using the relations: 



2tan(0/2) 
1 - tan 2 (fl/2) 

sn[F(0,s)] 
cn[F(0, s)j 



tan0 = 

tane = f!iSMJ, (15) 



Eq. (j!4p transforms into: 



sn(2w, s) _ 2 ' V «Xd 
cn(2r«7, s) -, _ /7 2 

V K^-cl 



w= ±y-u(r-r ). (16) 

The condition 4 ' \J\xdij) 7^ ±1 imposed by Eq. (fT6|) is equivalent to Xd{r)/a 7^ ±(2s 2 — 
l) -1 / 2 which is always fulfilled in the physical range of our interest. The l.h.s. in the first of 
Eq. (TTB1 can be rewritten using the double arguments relations: 



sn(2tu, s) 
cn(2zu, s) 



2sn(w, s)cn(zu, s)dn{w, s) 

1 — s 2 sn 4 {w 1 s) 
cn 2 (zj, s) — sn 2 (zu, s)dn 2 (zu, s) 
1 — s 2 sn 4 (zu, s) 



(17) 



and, through the definition in Eq. (Q, from Eqs. (fl6|) . (fT7|) I derive the two solutions for 
the classical equation of motion in the finite size model: 
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*2 } (r) 



4 AT" sn(va, s)dn(zj, s) 



2na- 



-W2kq 



cn[zu, s) 
cn(zu, s) 



sn(m, s)dn(zu, s) 



(18) 



Set now E = in Eq. (fTBj) . For k — 1/2, s 2 = 1 and the elliptic functions reduce to the 



limit cases 



3l|: 



sn(w, 1) = tanhw 
cn(zu, 1) = sechw 

dn(zu, 1) = sechw. (19) 



Then, from Eq. ( Tl8i) . I get consistently: 



s^(r) ±atanh[c;(r-r )/v / 2] 

-» Tacoth[cu(r - r )/v^]. (20) 

The first in Eq. (|20|) is the (anti)kink of Eq. ([8]) whose center of motion is set at x^(ro) = 
0. Instead, the second in Eq. fTSDl is the solution of Eq. (j7|) fulfilling the condition ^ (to) = 
oo. This proves the correctness of the analytical work. 

Computation of Eq. ([TBI) is an instructive task. Numerical investigation shows that 
both backgrounds predict the same physics. This is checked by establishing the one to 
one correspondence between E and the temperature T* at which the boundary conditions 
expressed by Eq. (J3J) are fulfilled: setting a value for E, both classical backgrounds fulfill 
Eq. (|5I) at the same T*. 

Note that both solutions present singularities. The zeros of (r), occurring at zu = 

(2) (2) 

±2nK(s) with integer n would be divergencies for x d (r). Viceversa, the zeros of x d (r) at 
w = ±(2n + l)K(s) are singular points for x d (t). 

For practical purposes, the choice of the first in Eq. (fT8]) turns out to be more convenient. 
In fact, taking the center of motion at x c /(r = 0) = 0, x d \r) properly describes a family 
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of finite size (anti)instantons which are bounded, continuous and odd in the range r G 
[-L/2,L/2\. 

The boundary conditions (see Eq. §5§) for x^\r) are always fulfilled before encountering 
the singular points, thus the paths are always continuous functions within their oscillation 



periods [l9j. Mathematically this means that, defining w* = ± 4 y/ r ^ujL/2, one finds numer- 



ically w*/K(s) < 1, hence the time period L is covered within a lattice period shorter than 
2K(s). 

Around x^\r), one can further study the finite size effects on the quantum fluctuation 
spectrum by means of the theory of the functional determinants. The details are given in 



Refs. 



321 ] and 33] while, in the next Section, I review the general features of the functional 



determinants method and evaluate the quantum fluctuations contribution for the case of the 
infinite size theory. 

III. Quantum Fluctuations around the (anti)kink 

For a particle with mass M moving in a potential V(x(t)) the probability amplitude 
to propagate, in the imaginary time L, from the initial position X{ (at —L/2) to the final 
position xj (at L/2) is given by the path integral [23] 



< Xf\Xi > L ~- 



.Dxexp 



A[x\ 



A[x] 



dr[-x{rf + V{x{r)) 

-L/2 \ 1 



(21) 



where A[x] is the Euclidean action that wejghs the contribution of the generic path x(r). 



Dx is the measure of the path integration 



341 ] . Accordingly, by expanding the Euclidean 



action in the neighborhood of the classical path, the propagation amplitude for a particle in 
the bistable quartic potential reads 
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< Xf\Xi >l= exp 

-L/2 



A[x d ] 
h 



Drjexp 



h 



A[x c 



-Zr/2 

L/2 



ilr[^xUr) + V(x d (r)) 



T 



r L/I { m i 

A/W= / dr ( Y V 2 (T) + -V"(x d (T)W 

V(x d (r)) = -^t{ x ci{t) -a) 2 (x d (r) + a) 2 
V"(x cl (r))=^(3x d (r)-a 2 ). 



Dr), defined below, is the integration measure over the path fluctuations. If L 
Eq. (jSJ) and the second of Eq. ( [22]) . one obtains the classical path action: 



(22) 
oo, from 



A[x 



cl\ 



3 6 • [Z6) 
This is also the (anti)kink energy in units of h exhibiting the peculiar dependence on 5 
emphasized in the Introduction. 

Taking the 5 — > limit (for a given uj) corresponds to decouple the two potential wells 
(a — > oo). Accordingly the classical energy required to connect the minima becomes infinitely 
large. A slightly different version of V{x d {r)) in Eq. (1221) is often found in the literature 
with uj replaced by In this case the action would be given by Eq. ( 1231) without the 

y/2 factor. 

The quadratic fluctuations contribution to the Euclidean action is derived by the third 
and fifth of Eq. ( !22l) . Performing an integration by parts on the kinetic term, the quantum 
fluctuation action transforms into: 



m r L ' 2 



dr 



-L/2 



d 2 r](r) 
dr 2 



+2cj 2 (1 - - cosh 2 vj)7](t)J7](t) 



Uj(r - T 



2. 



(24) 



The infinite size background of Eq. (jSJ) has been inserted in Eq. (!24j) . Then, consistently, 
Dirichlet boundary conditions T](t = ±L/2) = 0, with L — > oo, have to be assumed to 
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derive Eq. (pM)) from Eq. (|22|) . I'll come back to this point later. Expanding i](t) in a series 
of ortonormal components rj n (T) with coefficients <; n , 



(25) 



n=0 



one observes that the quantum fluctuation components obey a Schrodinger equation with 
Rosen-Morse potential 



Or] n (r) = e n r) n (r) 

d 2 3 
0= -— + 2u 2 (1 - - cosh" 2 w) 



(26) 



which admits two bound state eigenvalues (e , Ei) and a continuum spectrum (e n , n>2) 
J4]. In particular, the ground state 770 (t) is the Goldstone mode, discussed in the Introduc- 
tion, whose eigenvalue is £0 = 0. 

Thus, after solving Eq. (12"6"1) . the quantum action formally reads: 



M 



(27) 



n=0 



Note that, dimensionally: [e n ] = [sec 2 ], [r) n ] = [sec 1 / 2 ], [q n ] = [A sec 



-1/21 



Using the measure 



26|: 



/OO „ 
dv = n n / 
n=0 ~ 



fan 



(28) 



and performing Gaussian integrations over the directions q n one formally obtains the 
quantum fluctuations contribution to the propagator in Eq. ( |22l ): 



Drjexp 



h 



-1/2 



n=0 



(29) 
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The normalization constant K accounts for the Jacobian in the transformation to the 
normal mode expansion in Eq. (125]) . 

There is a singularity in Eq. (1291) caused by the e — eigenvalue which breaks the 
Gaussian approximation. This reflects the arbitrariness in the choice of r in Eq. (jSj). The 
latter represents in fact a one parameter family of solutions localized around their time 
center and, as such, they lack the time translational symmetry of the Hamiltonian. Any 
solution, in the set given by Eq. (jSj), can undergo an infinitesimal time translation which 
corresponds to an infinitesimal shift of To- The resulting shifted kink is however degenerate 
with the original one and therefore such a fluctuation costs no action. This is the physical 
origin of the zero mode. Accordingly we have the recipe to handle the divergence in Eq. (1291) : 
it can be easily seen that dx c i(r)/dT satisfies the Euler-Lagrange equation and solves the 
homogeneous differential problem associated to Eq. fl26|) . 0?7 n (r) = 0. Moreover it has no 
nodes being x c iij) monotonic. Then, dx c i(r)/dr is proportional to the Goldstone mode t} (t) 
of the ortonormal expansion in Eq. ( j25l) . 

These physical facts permit to tackle the divergence. In fact the integration measure in 
Eq. ( 128]) has a contribution along the % direction. Thus, transforming from the to the 



r -axis (with the norm of dx c i(r)/dT as Jacobian) one gets 26] 



7-co sjM/M V 2fi J V 21ft J_ Ln 



fL/2 

N- l =\l2l dr\x cl (r)\ 2 . (30) 



o 



The replacement (eq)^ 1 ^ 2 — > a/ 'A[x c i]/2irhL is often encountered in the literature. The 
latter however holds only in the L —>■ oo limit for which A[x c ;] = MN~ 2 whereas, for finite 
size/temperature systems, the replacement in Eq. ( 130]) is correct in general. 

Anyway the heart of the matter, unambigously shown by Eq. ( 130]) . is that the divergence 
in the path integral can be removed only by keeping L finite. It can be certainly large but 
it has to remain finite. On the other hand the (anti)kink solution in Eq. (jSj) is intrinsically 
based on the assumption L — >• oo. Moreover, as emphasized above, its derivative (oc ffafr)) 
satisfies Dirichlet boundary conditions only at L — > oo. One may face the contradiction 
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by forcing the fluctuation to vanish in a finite r— range: this would compress the particle 
distribution thus lifting the Eq eigenvalue upwards and removing the zero Although the 
discrepancy (with respect to the asymptotic picture) is of order 0(exp(— uL)) and therefore 
it may be controlled at large L, I felt that there was a need for a consistent formulation of 
the theory able to incorporate finite size effects already at the classical stage. This justifies 
the work presented in the previous Section. 

Hereafter I complete the derivation of the path integral based on the infinite size 
(anti)kink. Thus, Eq. (|29|) can be rewritten as: 



Dt]{t) exp 
F{uL) 



My) 
h 



f \ V 2 

F(uL) {i.L=o £ 7 



n 



n=2 £ n 



1/2 
cont 



I ll n =0 fc r 



1/2" 



(31) 



Formally, I have multiplied and divided by the product of the harmonic quantum fluc- 
tuation eigenvalues [251 ] . F(luL) is the well known harmonic fluctuation factor which 
absorbs K 25] . The two bound states have been factored out so that the set of continuum 
eigenvalues (denominator) is discrete but it contains two states less than the harmonic spec- 
trum (numerator). Then, the evaluation of the quadratic quantum fluctuation contribution 
is practically reduced to the estimate of the ratio of functional determinants in Eq. (I3T1) . 



This has been done by several methods in the literature 
ideas by Gelfand and Yaglom 



35 



361 . 137] extending fundamental 



. The fact that a ratio of determinants 



34j and Forman 

appears is not a simple mathematical trick. In fact each determinant, being the product 
over an infinite number of eigenvalues with magnitude greater than one, would diverge in 
the T — > limit. Such (exponential) divergences cancel out in the ratio which is meaningful 
both in value and sign and arises naturally in the path integral method |34|. 

In the instantonic approach, the essential observation is that both determinants satisfy the 
Jacobi equations (that is, the homogeneous differential equations associated to the respective 
stability equations) around their classical solutions, see Ref. 26| for an extensive discussion. 
The remarkable path integral result which follows is that the ratio is determined only on 
the base of the classical action (or of the classical path velocity) and it yields a factor 9a; 4 . 
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Then, including also the bound state eigenvalue e\ = 3u; 2 /4 in Eq. (I3~TT) and taking the 
square root of the whole, the fluctuations contribution gets: 



Di](t) exp 



MN~ 2 
2V3cuF(tuL)\l^- f -L. 



(32) 



This result completes the calculation of the quantum fluctuation for the infinite size case 
in the semiclassical approximation of Eq. (122]) • 

IV. Infinite and Finite Size Path Integrals 

Assume now that, in the propagation amplitude, the initial and final positions are very 
close to the potential minima: Xi ~ —a, Xf ~ a and say ipo(x±a) the ground state harmonic 
wavefunctions in the potential wells. For x ~ ±a, (po(x+a)(po(x— a) ~ \J ^r. The harmonic 
factor, F(luL) can be approximated in the large L limit as, F(luL) ~ yl exp(—ioL/2). 
Then, from Eq. ( 122]) and Eq. ( |35l) . the path propagator for a particle in a single double well 
potential reads: 



< -''./ 1 >l- n = \/ exp (-uL/2)tt 00 L 



{loo = exp 



A[x d ] 

h 



9 R 1 MN ~ 2 
2 V Sou 



2nh 

(33) 



This is the general result of the infinite size instantonic theory. The quantum tunneling 
removes the two fold ground state degeneracy and induces a splitting in the energy levels. 
The size of this quantum effect is governed by the classical action hosted in the expression of 
Qoo. The classical path solution determines the shape of the potential in the stability equa- 
tion (Eq. (12"6T) ) thus providing the frame within which the quadratic quantum fluctuations 
display their strength and set the tunneling rate. 

Finally and fundamentally, note that does not depend on temperature in the infinite 
size instantonic approach as A[x c /] is a constant in such limit. This limitation of the theory 
is removed in the finite size method based on the classical background of Section II. 



14 



Around such background one has to evaluate the quantum fluctuations contribution by 



finding the functional determinants ratio for a finite size system [35j after extracting the zero 
mode as shown above. Working out the calculations, I get the semiclassical path integral 
for one (anti)instanton in finite size L: 



< Xf\Xi > L 

Q{L) = 

n QF (L) 



M 



2nhL 2sinh(o;L/2) 



n(L)L 



cxp 



A[x d ] 

h 

MN-- 
2ixh 



n QF (L) 

Det[h] 



Det R [0} 



(34) 



with: 



Det[h] = -4sinh 2 (a;L/2) 
h = -d 2 T +u 2 

oo 

Det R [6) = l[e n (35) 

n= 1 

Note that the harmonic determinant has been, in turn, normalized over the free particle 
determinant Det[—d 2 ] = J which incorporates the constant K. In the L — >• oo limit, I 

get 

Det[h ] -+ 12u* (36) 
Det R [0] 

thus recovering the value of the infinite size instantonic approach given in Eq. (1551) . This 
proves the correctness of the analytical procedure. 

Q® F (L) in Eq. (134]) accounts for the quantum fluctuation effects in the finite size tunneling 
problem while Q(L) is the overall tunneling frequency which removes the twofold degeneracy 
of the double well potential in a finite domain. From the path integral in Eq. (1341 one can 
extract the physical properties of the closed quantum system for any size. 
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V. Final Remarks 



I have reviewed the instantonic method for the propagation amplitude of a particle moving 
via tunneling in a bistable quartic potential. This has been done within the standard 
semiclassical scheme for a 4 model admitting the (anti)kink solution as classical path which 
interpolates between the vacua in a infinite time. Since the time is conceived as imaginary by 
a Wick rotation, an infinite time description maps onto the zero temperature limit according 
to the Matsubara formalism. The classical action and the stability equation for the quantum 
fluctuations have also been solved in the asymptotic limit. However, in deriving the particle 
propagator I have pointed out that, at some steps, the time (the size of the system) is in 
fact assumed to be large but finite in order to tackle the divergency due to the quantum 
fluctuation zero eigenvalue related to the translational invariance of the system. Led by this 
observation, I have described the general method which permits a consistent extension of the 
instantonic approach such to incorporate finite size effects both in the classical background 
and in the quantum spectrum. The fundamental step, the replacement of the infinite size 
(anti)kink by the finite size classical path of Eq. ffTSj) . has been explained in detail making 
use of the formalism of the Jacobi elliptic functions. Along similar patterns one generalizes 
the bounce solution of the Euler-Lagrange equation for metastable potentials to the finite 
size case [l]]. As a general consequence of the finite size formalism, the quantum fluctuation 
spectrum obtained in semiclassical approximation becomes dependent on the system size or, 
in the Matsubara formalism, on the system temperature. From a physical point of view, 
the elliptic functions formalism permits to evaluate the size/temperature effects both on the 
tunneling energy of bistable potentials and on the decay rate in systems with metastable 
potentials. 
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